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NUMERICAL  AND  ANALYTICAL  STUDIES  OF  BEAM  CHANNEL  TRACKING 


I.  Introduction 

An  intense  relativistic  charged  particle  beam  may  be  injected  into  a  pre¬ 
formed  channel  in  a  neutral  gas.  The  channel  may  have  been  formed  in  a 
variety  of  ways,  i.e.,  with  a  laser,  an  electrical  discharge,  or  a  previous 
particle  beam,  and  it  may  be  simply  a  heated  reduced-density  region  in 
pressure  balance  with  the  surrounding  neutral  gas  (a  density  channel),  or  it 
may  be  ionized  (a  conductivity  channel),  or  both.  If  the  beam  is  displaced 
from  the  axis  of  the  channel,  the  question  arises  as  to  whether  electro¬ 
magnetic  forces  are  generated  which  tend  to  restore  the  beam  to  the  channel 
axis  or  push  it  further  away.  We  shall  refer  to  these  as  channel  tracking  or 
detracking  forces,  respectively. 

.In  this  paper,  we  report  on  theoretical  calculations  of  channel  tracking 
phenomena  for  the  case  where  there  is  no  pre-existing  current  flowing  in  the 
channel.  In  this  situation,  the  question  of  whether  the  beam  will  stay  in  the 
channel  separates  into  three  distinct  problems.  The  first  of  these  relates  to 
the  forces  between  the  beam  and  the  channel.  If  the  pre-existing  channel 
conductivity  a(r)  is  not  too  high  (4itar^  <  c,  where  r^  is  the  beam  radius  and 
c  is  the  speed  of  light),  then  electrostatic  forces  are  slightly  dominant  in 
the  front  of  the  beam,  back  to  about  the  pinch  point^.  Since  electrostatic 
image  forces  are  predominately  attractive,  one  might  expect  the  front  of  the 
beam  to  be  attracted  toward  the  channel.  However,  the  actual  situation  turns 


complicate  matters.  We  shall  refer  to  these  questions  as  the  tracking 
problem,  proper.  This  alone  is  the  subject  of  this  paper. 

In  the  main  part  of  the  beam,  well  behind  the  pinch  point,  the  gas 
conductivity  is  large  enough  to  insure  space  charge  neutrality,  and  thus 
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magnetic  forces  predominate.  The  force  between  the  beam  and  the  oppositely- 
directed  plasma  current  in  the  channel  is  repulsive,  and  generally  much 
stronger  than  the  forces  on  the  beam  head.  Nonetheless,  the  beam  is  not 
necessarily  pushed  out  of  the  channel,  because  there  are  magnetic  coupling 
forces  that  attract  the  beam  toward  the  beam  axis  established  by  the  head, 
which  is  itself  attracted  toward  the  channel  (under  some  circumstances).  The 
question  of  whether  these  beam-beam  forces  are  strong  enough  to  balance  the 
repulsive  beam-channel  forces  and  permit  an  equilibrium  to  exist  will  be 
referred  to  as  the  self-coupling  problem.  The  equilibrium,  if  it  exists,  may 
be  such  that  the  displacement  of  the  beam  axis  from  the  channel  axis  increases 
with  5  S  ct  -  z,  the  distance  behind  the  beam  head.  This  will  be  considered 
in  a  subsequent  paper. 

Assuming  that  an  equilibrium  exists,  the  question  remains  as  to  its 

stability  properties,  particularly  to  the  hose  instability.  This  stability 

2 

problem  has  been  studied  at  great  length  ,  but  only  in  regard  to  a  beam  in  a 
uniform  medium,  or  a  beam  whose  equilibrium  is  assumed  to  be  on-axis  in  a 
channel.  It  is  well-known  that  the  presence  of  plasma  return  current  is  a 
destabilizing  factor.  However,  stability  properties  remain  a  subject  for 
future  study,  as  regards  the  correct  beam  equilibria  in  the  presence  of  a 
density  or  conductivity  channel. 

In  this  study,  we  consider  a  wide  range  of  possible  channel  states.  For 


specificity, 


the  gas  is  taken  to  be  air  at  1  atm  outside  the  channel,  with  on- 
density  ranging  from  0.1  atm  to  1  atm,  and  a  simple  model  for 
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ionization  by  the  beam  is  included  in  the  calculation.  Within  our  model,  the 


channel  density  affects  the  forces  on  the  beam  only  indirectly,  by  means  of 
the  density-dependence  of  avalanche  ionization  and  recombination,  which  in 
turn  affect  the  conductivity  of  the  gas.  The  pre-existing  channel 
conductivity  is  varied  from  zero  to  a  large  value.  The  channel  density  is 
assumed  to  be  independent  of  time  and  of  the  axial  position  z;  the  initial 
conductivity  is  also  independent  of  z,  but  o  does  depend  on  ?  due  to  beam- 
produced  ionization.  Our  model  also  assumes  that  the  initial  o(r,5)  is 
axisymmetric  about  the  original  channel  axis.  The  radial  profiles  of  the 
initial  channel  density  and  conductivity  are  assumed  to  be  Bennett. 

In  a  real  experiment,  the" ratio  of  the  beam  radius  r^  to  channel  radius 
rc  varies  quite  a  lot  from  the  head  of  the  beam  to  the  tail.  Near  the  head, 
the  beam  has  a  trumpet  shape*,  so  its  radius  may  be  much  larger  than  the 
channel  radius.  Near  the  tail,  the  beam  is  well  pinched  so  its  radius  can  be 
smaller  than  the  channel  radius.  Ideally  one  should  use  a  self-consistent 
envelope  equation  for  the  beam  radius  but  we  have  not  incorporated  that  in  our 
model  at  this  point.  We  have  adopted  a  simpler  approach,  which  is  to  use  a 
uniform  pencil  beam  and  study  the  dependence  of  tracking  force  on  the  ratio  of 
r^/rc.  Occasionally,  we  do  assign  an  envelope  shape  to  the  beam,  but  this 
envelope  is  not  coupled  back  to  the  fields;  therefore,  we  should  look  at  those 
results  as  only  qualitative. 

Throughout  each  run  the  beam  is  modeled  as  a  rigid  straight  rod  which  may 
be  displaced  from  the  channel  axis  by  a  uniform  time-independent  displacement 
Y.  We  emphasize  that  we  do  not  assume  that  Y  is  small.  The  geometry  is  shown 
in  Fig.  1.  Because  neither  the  beam  nor  the  channel  are  allowed  to  respond 
dynamically,  we  can  calculate  unambiguously  the  tracking  force  between  the 


resultant  force  acting  on  the  beam  would  have  an  additional  component  due  to 
the  coupling  force  between  beam  slices.  This  force  is  typically  at  least  two 
orders  of  magnitude  larger  than  the  electrostatic  tracking  force  and  therefore 
the  delicate  physics  of  beam  channel  tracking  can  easily  be  overwhelmed.  The 
adoption  of  a  time  independent  approach  also  allows  us  to  eliminate  the  hose 
motion  (instability)  which  does  not  concern  us  in  this  report. 

The  numerical  study  is  performed  using  the  DYNASTY  II  code  developed  at 
NRL.  It  solves  Lee's  ultrarelativistic  approximation  to  the  Maxwell  field 
equations^, 


where  Az  is  the  axial  component  of  the  vector  potential,  $  is  the  scaler 
potential,  and  a  is  the  gas  conductivity,  together  with  an  air  conductivity 
model  that  includes  direct  ionization,  avalanche  and  recombination, 

tf(r,S)  -  <Jb(r)  +  (aE/c)o  -  Sra2.  (2) 

-4  -15 

Here  <  m  8.8x10  cm/statcoul,  8  ■  7.1x10  (p  /p  )sec/cm,  the  avalanche 

r  go 

coefficient  av  is  a  function  of  the  ratio  of  the  electric  field  E  to  the  gas 

density  p  ,  given  as  Eqs.  (17)  and  (18)  of  Ref.  5,  and  p  is  the  air  density 
g  o 

at  1  atm.  DYNASTY  II  is  not  linearized  about  cylindrical  symmetry,  as  are 
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typical  monopole/dipole  beam  codes^’^;  therefore,  we  are  not  restricted  to 

study  only  small  beam  displacement  as  required  by  linearized  codes.  This 

capability  has  allowed  us  to  discover  a  new  beam-channel  equilibrium  which  was 

not  observed  before.  To  elucidate  our  rather  surprising  numerical  results,  we 

have  also  analyzed  a  solvable  two-dimensional  electrostatic  tracking  model 

which  yields  similar  conclusions  about  the  beam  equilibria  with  respect  to  the 

channel.  The  stability  and  propagation  characteristics  with  respect  to  these 

equilibria  are  unknown  at  this  point. 

To  summarize  our  main  results,  we  fiv.'d  that  the  beam  head  always  tracks 

the  channel  provided  the  preformed  channel  conductivity  a  (■  4irar, /c)  is  not 

o  b 

bigger  than  0.8  and  not  much  smaller  than  0.1.  For  a  pencil  beam,  the 

equilibrium  position  of  the  beam  head  can  be  on  or  off  the  channel  axis 

depending  on  the  ratio  of  the  width  of  the  beam  and  the  channel.  A  pure 

density  channel  results  in  a  beam  equilibrium  off-axis  in  the  channel  for  all 

1  • 

cases  studied.  The  body  of  the  beam  behind  the  pinch  point  always  feels  a 
detracking  force,  and  therefore  will  stay  in  the  channel  only  if  the  self¬ 
coupling  force  holds  it  to  the  head.  The  electrostatic  tracking  force  is  very 
weak —  as  much  as  five  orders  of  magnitude  weaker  than  the  pinch  force  at  full 
current,  for  typical  cases  studied  with  rise-time  for  the  current  in  the 
range  30  r^/c  to  60  r^/c.  For  typical  laboratory  experiments^’^,  this  force 
is  too  small  to  be  detected. 

Section  II  discusses  the  numerical  results;  Section  III  describes  the 
analytical  calculations;  and  Section  IV  presents  the  conclusions. 


II.  Numerical  Results 

In  this  section,  we  present  numerical  results  for  the  dependence  of  the 
tracking  force  on  the  beam  radius,  channel  radius,  channel  pressure,  channel 
conductivity,  and  current  rise  time.  The  geometry  of  the  system  is  shown  in 
Fig.  1  where  the  preformed  channel  is  displaced  as  a  rigid  rod  from  the  beam 

a 

axis  by  an  amount  e  .  The  force  on  a  slice  of  the  beam  at  £  is  given  by 

"y 

Ft(0  -  Kq  J  d2r  Jb(e)  [Az(0  -  $(?)]  *  S  »  W 

y 

where  is  the  beam  current  density.  The  beam  and  the  channel  profiles  are 
chosen  to  be  Bennett.  Kq  is  a  normalization  constant,  chosen  to  normalize  the 
tracking  force  to  the  pinch  force  Fp  of  a  beam  at  full  current.  We  fix  the 
beam  radius  to  be  r^  *  0.5  cm  for  pencil  beams.  The  channel  radius  is  varied 
from  0.5^  to  1.5r^.  All  lengths  (?  and  rc)  are  normalize^  to  r^,  and 
conductivity  is  expressed  in  the  dimensionless  units  4irorb/c.  The  maximum 
beam  current  is  I^Q  *  10  kA  and  it  ramps  up  in  5  according  to 

Ib  -  IbQ  [1  -  exp  -U/4r)2],  (5) 

where  ;  =  ct  -  z  is  a  measure  of  the  distance  along  the  axis  from  the  beam 
head. 

We  have  checked  out  the  numerical  validity  of  the  code  in  the  following 
ways.  For  Y  «  1,  we  compare  the  results  with  the  NRL  linearized  code  SIMMl^, 
and  find  the  forces  calculated  by  the  two  codes  differ  by  only  a  few  percent. 
For  Y  >  1,  we  alternately  displace  the  beam  and  the  channel  by  the  same  Y  to 
make  sure  they  give  the  same  result  (within  a  few  percent  error).  As  an 
additional  check,  we  displace  both  the  beam  and  the  channel  together  by  the 
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same  Y  and  we  find  the  force  (which  should  be  zero)  is  typically  less  than  one 
percent  compared  to  when  just  the  channel  is  displaced. 

We  shall  use  the  following  abbreviations  to  label  the  channels: 
a  channel  =  a  pure  conductivity  channel, 

n  channel  =  a  pure  density  channel, 

a  +  n  channel  =  a  combination  of  conductivity  and 

density  channel. 

The  first  set  of  runs  are  for  pencil  beams,  i.e.,  r^  constant. 

Typical  ^-dependences  of  beam  current  density  on  the  beam  axis,  plasma 
current  density  J_  on  the  channel  axis,  axial  electric  field  E_  on  the  channel 
axis,  and  conductivity  a  on  the  channel  axis  near  the  beam  head  are  shown  in 
Fig.  2.  Figures  3-6  show  the  dependence  of  the  tracking  force  Ft(s)» 
normalized  to  pinch  force  Fp  at  full  current,  on  the  channel  radius  rc  and 
displacement  Y.  Figure  3  shows  plots  for  rc  *  2.0  and  Y  «•  0.1.  The  a  and 
a  +  n  channels  show  tracking  while  the  n  channel  does  not.  The  peak  tracking 
force  for  the  a  and  a  +  n  channels  are  about  the  same.  These  findings  are  in 
agreement  with  previous  calculations  that  were  linearized  by  assuming  Y  <<  r^ 
and  Y  «  rc^®.  The  a  channel  initially  has  oq  »  0.1  on  axis  and  the  n  channel 
has  density  0.1  atm  on  axis.  If  the  displacement  is  increased  to  Y  *  1.5, 
while  keeping  rc  »  2.0  (Fig.  4),  the  n  channel  now  shows  relatively  good 
tracking  (about  ten  times  stronger  than  the  a  channel  with  Y  ■  0.1).  The 
tracking  force  for  the  a  and  a  +  n  channels  is  even  stronger,  by  another  order 
of  magnitude.  Figures  3  and  4  show  that  a  and  a  +  n  channels  allow  the  beam 
to  track  on  axis,  with  the  tracking  force  stronger  for  a  a  +  n  channel  than 
for  any  other  channel  type,  while  the  beam  in  an  n  channel  feels  a  tracking 
force  only  when  it  is  off  axis.  This  off-axis  tracking  implies  that  the  beam 
head  is  in  unstable  equilibrium  when  it  is  propagating  slightly  off  axis  in  a 


channel,  but  that  an  equilibrium  exists  with  the  beam  head  off-axis  by  a 
prescribed  amount.  Moreover,  the  equilibrium  value  of  Y  is  an  increasing 
function  of  It  should  be  emphasized  that  these  results  cannot  be  obtained 
with  linearized  codes  because  these  codes  are  restricted  to  values  of  Y  <<  1. 

Figure  5  shows  similar  tracking  force  plots  for  r£  =  0.5.  These  results 
should  give  us  an  indication  of  how  well  the  beam  tracks  near  the  blown-up 
head  region.  We  find  that  the  beam  feels  a  detracking  force  in  any  type  of 
channel  for  values  of  Y  up  to  1.0  (curves  a-c).  For  larger  displacement,  Y  * 
1.5,  the  force  is  reversed  and  tracking  occurs  (curves  d-f).  These  results 
suggest  that  the  blown-up  beam  head  will  never  propagate  coaxially  with  the 
preformed  channels,  but  can  reach  an  equilibrium  position  off  axis.  Whether 
this  will  drive  hose  instability  remains  to  be  determined. 

In  Fig.  6,  we  show  the  dependence  of  the  tracking  force  as  a  function  of 
displacement  Y  at  a  fixed  value  of  £  -  9.0,  for  a  case  with  rc  *  0.5rb.  This 
clearly  Indicates  that  the  beam  is  unstable  to  small  displacement,  but  that  at 
larger  Y,  the  beam  polarizes  the  channel  and  the  force  becomes  attractive. 
Consequently,  for  rc/rb  <  1.0,  the  beam  head  will  sit  somewhere  off  axis  with 
respect  to  the  channel. 

The  cases  for  r£  *  1.0  are  shown  in  Fig.  7.  This  situation  resembles  the 
neck  region  of  the  beam.  Again  we  find  for  all  types  of  channels  that  the  beam 
detracks  for  small  Y  and  tracks  for  larger  Y,  with  the  a  +  n  channel  giving 
the  strongest  tracking.  Notice  that  in  contrast  to  the  case  with  r£  »  0.5 
(Fig.  5),  the  case  Y  ■  1.0  now  gives  some  tracking.  This  indicates  the  larger 
the  channel  radius,  the  easier  it  is  for  the  beam  to  track  the  channel.  This 
point  is  illustrated  in  Fig.  8  where  we  compare  the  tracking  force  for  cases 
with  rc  *  2.0  and  5.0,  with  Y  ■  0.1  in  both  cases.  The  case  with  r£  =  5.0 


V*. 


gives  a  much  stronger  tracking  force  and  its  tracking  region  in  ?  is  twice  as 
long. 

The  dependence  of  the  tracking  force  on  the  current  rise  time  X,  ^ 

(measured  in  units  of  £)  is  shown  in  Fig.  9.  A  longer  current  rise  time 
essentially  "stretches"  the  tracking  region  in  5  and  weakens  the  peak  tracking 
force,  but  it  does  not  change  the  sign  of  the  tracking  force. 

For  all  the  cases  presented  in  Fig.  3-9,  beam-channel  forces  become 
detracking  at  larger  £.  This  is  to  be  expected  because  in  the  magnetic  regime 
(where  the  gas  conductivity  is  high  and  the  beam  space  charge  is  neutralized), 
the  repulsive  magnetic  force  between  the  beam  and  the  return  current 
dominates.  This  does  not  mean  that  part  of  the  beam  is  repelled  and  breaks 
away  from  the  beam  head.  Once  the  beam  head  tracks  the  channel,  the  self- 
coupling  force  is  usually  strong  enough  to  hold  the  rest  of  the  beam  together 
provided  the  initial  conductivity  is  not  too  large.  The  dynamical  response  of 
a  beam  to  various  types  of  channels  in  the  linearized  regime  has  been  studied 
by  Masamitsu,  et  al.^®. 

There  exists  a  critical  conductivity  oq  ■  c?c  above  which  the  beam  never 

tracks  the  channels  for  any  values  of  r.  and  Y.  For  a  >  a  ,  the 

c  o  c 

electrostatic  regime  is  compressed  to  such  a  small  region  as  to  render  channel 

tracking  totally  ineffective.  In  the  other  extreme,  if  oq  «  0.1  and  there  is 

no  density  channel,  the  tracking  force  becomes  very  weak  and  linearly 

proportional  to  oq  for  the  cases  of  on-axis  equilibrium^.  For  cases  we  have 

studied,  is  around  0.8,  and  the  optimum  conductivity  for  on-axis  tracking 

is  around  a  m  0.1. 
o 

We  have  not  yet  included  a  self-consistent  envelope  equation  in  our  code, 
but  we  can  estimate  the  effect  of  a  trumpet-shaped  beam  head  on  tracking  by 
assigning  an  envelope  to  the  beam  radius  (normalized  to  the  fully-pinched  beam 
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•N 


radius) , 


rb  -1.0+8  e‘?/ct,  (6) 

where  we  choose  8  -  10.0  and  a  -  2.0.  We  fix  r£  -  1.0.  Consistent  with  the 
narrow-channel  pencil  beam  case  (Fig.  5),  the  beam  is  found  to  detrack  for  Y 
«  1.0  and  track  for  larger  displacements.  Figure  10  shows  the  tracking  force 
in  <j,  a  +  n  and  n  channels  with  Y  -  1.5.  The  strongest  force  occurs  in  a  +  n 
channel.  Channel  tracking  of  a  beam  with  a  self  consistent  envelope  will  be 
reported  in  a  later  publication. 

Figure  11  shows  the  variation  of  the  n  channel  force  for  two  cases  with 
Y  -  1.5  and  on-axis  channel  density  p  (r  -  0)  »  0.1  atm  and  0.2  atm.  In  each 

s 

case,  the  density  outside  the  channel  is  1  atm.  This  clearly  indicates  a 

deeper  density  channel  gives  better  off-axis  tracking,  as  might  be  expected. 

« 

The  results  for  the  parametric  dependence  on  r£,  Y  and  types  of  channel 
by  the  tracking  force  are  summarized  in  Table  I  for  the  pencil  beam  and  Table 
II  for  the  trumpet  beam.  The  last  column  represents  the  maximum  tracking 
force  normalized  to  the  pinch  force. 

In  Fig.  12,  we  show  the  resultant  tracking  force  ft  on  a  ring  of  beam 
electrons  located  between  r  and  r  +  Ar  at  a  given  c,  for  Y  -  0.1  and  1.5. 

A  a  +  n  channel  is  chosen  for  this  example.  For  both  values  of  Y,  the 
central  core  of  the  beam  always  feels  a  net  force  that  attracts  it  toward  the 
channel,  even  though  the  net  force  on  the  beam  slice  as  a  whole  is  detracking 
for  Y  -  0.1.  This  phenomena  has  previously  been  observed  in  linearized 
simulation  codes  using  wide  channels^.  However,  the  portion  of  the  beam  that 
feels  an  attractive  force  is  very  small  (r/r^  <  1/3)  for  the  Y  -  0.1  case. 
Since  the  tracking  or  detracking  force  is  usually  10  to  10  J  times  smaller 


kV*T 


III.  Analytic  and  Model  Calculations 

The  DYNASTY-II  code  calculations  presented  in  previous  sections  have 
arrived  at  conclusions  regarding  the  beam-channel  interaction  that  were  not 
anticipated  and  may  seem  peculiar,  namely  that  if  the  channel  is  initially 
weakly  ionized  (4ira  r^/c  <  0.8)  the  beam  head  is  attracted  toward  the  channel 
axis  if  the  beam  is  narrow,  but  is  repelled  from  the  channel  axis  and 
attracted  toward  a  certain  fixed  displacement  Y  *  YQ  if  the  beam  is  broad. 

(If  the  channel  is  initially  unionized,  the  beam  is  practically  force-free 
for  Y  <  Yq,  but  feels  a  tracking  force  when  Y  >  YQ,  whether  r|j/rc  is  large  or 
small.)  In  this  section  we  shall  consider  a  simple  example  which  can  be 
explored  by  means  of  physical  arguments  and  elementary  mathematics,  and  which 
elucidates  the  origin  of  these  effects. 

We  consider  a  beam  of  specified  charge  density  profile  p,(?,  r)  and  we 

D 

consider  only  the  limit  of  small  o,  applicable  to  the  head  of  a  beam 
propagating  in  neutral  gas  or  in  a  weakly  pre-ionized  channel.  It  is  worth 
mentioning  that  for  an  ultrarelativistic  beam  (y  »  1)  Maxwell's  equations 
reduce  exactly  to  Lee's  field  equation  (1)  in  the  limit  of  small  o  (as  well  as 
the  limit  of  large  o  and  several  other  limits)^.  It  has  further  been  pointed 
out  that  to  first  order  in  o  the  forces  on  the  beam  due  to  Eq.  (1)  reduce  to 
two  dimensional  electrostatics^.  We  shall  reproduce  the  expansion  that  leads 
to  this  conclusion. 

To  order  o°,  Eq.  (lb)  yields 


7*  A(o)  -  V*  *(o)  +  *(Ei), 


where  tp  is  some  function  independent  of  £.  Since  both  A  and  <j>  are  zero  ahead 
of  the  beam  (no  precursor  signals  can  get  ahead  of  the  beam  in  the 
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ultrarelativistic  limit),  and  A  and  $  can  be  taken  to  be  equal  at  large  radii 
(either  metallic  or  vacuum  boundary  conditions  at  large  r),  we  find 


A<°>  -*(0), 


From  Eq.  (la) 


-  4*  PV 


i  •  g  •  y 


*(0)(r)  -  2  /  —  /  dr  2-rrr  p  (r  ), 


r  r  o 


where  is  the  beam  charge  density.  To  order  a°,  the  bear  is  force-free: 
its  vacuum  electric  and  magnetic  self-forces  cancel  (in  the  ultrarelativistic 


limit). 


To  next  order,  Eqs.  (la)  and  (8)  with  appropriate  boundary  conditions 


yield. 


L<1>  - 


and  Eq.  (lb)  yields, 


V2  ,  -  A  ,  „(D 

s>  ■  -  4  ir  p  , 


where  p  is  the  plasma  charge  density,  given  by, 


i  •  *  •  *  *  *  «  *  •  *  '  i  *  «  "  «  "  •  •  ■  •  •  •  *  •  *  >.VL  1  «  •  "  •  1  •  'V  *•  *  *  •  .  ,  ,  ’ 


m 


Through  first  order  in  a,  the  force  on  the  beam  is  entirely  due  to  the 
electrostatic  potential  <j>^,  which  is  specified  by  the  two-dimensional 
electrostatic  field  equation  (12),  with  the  plasma  charge  distribution  (13)  as 
the  only  source.  Henceforth,  we  simplify  notation  by  dropping  the 
superscripts  on  and  p^\  Equations  (7)-(13)  have  previously  been 

derived  by  Briggs  and  Lee**. 

Next  we  consider  the  specific  case  of  a  pencil  beam  of  Gaussian  profile. 


pbo  eXp(_  r2/rS)» 


propagating  in  a  specified  displaced  Gaussian  conductivity  profile. 


a  -  <JQ(? )exp[- — +  Y)-"]  =  0Q(c)exp(-  ^j-). 

rb/a  rb 


Equations  (10)  and  (13)  may  be  integrated  in  closed  form  for  this  case, 
yielding 


2 

16ir  p.  ,  ,r  2 ,  -2 

— ~  /o  d?  ff0(?  ) ft 1  +  °  +  5  exp( " 

r  rb 


«(*  +  exp(" 

r  r.  J 


We  could  solve  Poisson's  equation  (12)  and  calculate  the  resultant  force  on 
the  beam  for  any  Y,  but  this  could  only  be  carried  out  numerically  and  would 
not  be  particularly  illuminating.  We  therefore  follow  the  alternative  course 


\.*v\vv  ^ 


*  .  •  ■  a  m 


*•••••>  •  **  • 
--••.■•..-•.s'. 


,*•  y*  *"•  ^ 


,  •  a  •  .  !  a  f  -  a  m  m  • 
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of  examining  separately  the  limits  of  small  Y  and  large  Y. 


If  Y  «  r^  and  we  linearize  in  Y,  Eq.  (16)  may  be  written  as  the  sum  of  a 
term  symmetric  about  the  conductivity  channel  axis  and  another  term  symmetric 
about  the  beam  axis, 


p(Si)  -  p1(r;a)  +  P2(r;a), 


(17) 


where  r  is  the  distance  from  channel  axis, 


Px(r;a)  -  -  a  16ir  pfeo 


-  {  ~2 
/q  d?  a(0  <  “  exp(- 


+  exp 


-2 

(1  +  a)r 

4 


*  T  El  -  T  E1 

rb 


(1  +  a)r‘ 


-2  -i 


(18) 


and 


2  eC  *■ 

,(r;a)  -  -  16ir  pbo  d?  o(?)exp 


[^1 


Ei(“f)  +f EiTii-Lfoi2]  , 

rL  L  rL  -1 


(19) 


where  Ei(x)=  -  /  dtt  e  is  the  exponential  integral.  Because  of 
-x 

axisymmetry,  fields  derived  from  p2(r)  exert  no  resultant  force  on  the  beam, 
and  thus  the  resultant  force  on  any  given  beam  slice,  to  first  order  in  Y  (or 
the  total  potential  of  a  beam  slice  to  second  order  in  Y)  derives  entirely 

v  v 

from  p^(r).  Since  p^  is  itself  axisymmetric  as  a  function  of  r,  the  potential 
associated  with  the  charge  distribution  p^  is  immediately  specified  to  be 


”/  **tt*  r»#. 


<t> L( r;  a) 


rr  dr  rr  .  ", 

“  J o  —  /0  dr  4irr  pl(r  }* 
r 


We  note  from  Eq.  (18)  the  following  properties  of  p. 


PbQ  Px(r)  <0  for  0  <  r  <  rQ, 


(21a) 


Pbo  PL(r>  >  o 


for  r  <  r  <  •  , 
o 


(21b) 


Pbo  /  dr  2irr  p^r)  *  4irpbo  /  d?  p(?)  >  0, 


(21c) 


where  rQ/rb  is  a  function  of  a  determined  by  the  equation, 


exP(-  ar^/rj)  -  I  Ei(a  r*/r£)  -  e*p[-  (1  +  a)r*/r*]  -  j  Ei[(l  +  a)r*/r£]. 


(2  Id) 


Properties  (21a)  and  (21b)  might  well  be  expected,  since  (for  an  electron 
beam,  with  pbQ  <  0)  a  positive  image  charge  collects  in  the  region  of  highest 
conductivity,  leaving  a  negative  plasma  charge  elsewhere.  Equation  (21c)  may 
be  a  bit  surprising,  since  the  plasma  must  be  globally  neutral;  however  it  is 
another  manifestation  of  image  charge  formation:  the  component  p^  has 
negative  total  charge,  while  the  component  p^,  which  is  closer  to  the  beam  but 
which  exerts  no  resultant  force  on  the  beam,  has  positive  total  charge.  These 
properties  may  reasonably  be  expected  to  hold  generally  for  a  beam  with  any 
monotonic  profile  p,  (r)  in  a  channel  o(r)  of  monotonic  profile. 

D 

We  consider  the  forces  on  a  beam  slice  displaced  from  the  axis  of  the 
central  potential  <|» ^ ( r )  of  Eq.  (20).  We  shall  consider  the  beam  slice  to  be  a 


rigid  object,  made  up  by  fixed  point  elements  of  charge  density  or 

of  fixed  rings  of  charge  density  at  radius  r.  The  force  on  a  ring  is  the 
resultant  of  the  forces  on  its  elements;  the  force  on  the  whole  beam  slice  is 
the  resultant  of  the  forces  on  its  rings.  No  account  is  taken  here  of  the 
fact  that  the  beam  is  made  up  of  moving  electrons  that  are  not  confined  to  any 
particular  ring,  or  that  the  beam  may  tend  to  distort  from  its  circular  shape, 
or  even  tear  apart,  because  the  channel  force  is  different  on  different  parts 
of  the  beam.  In  other  words,  we  are  only  calculating  the  forces  on  the  beam 
here,  not  following  the  beam  dynamics  at  all. 

Using  the  channel-centered  coordinate  system  r,  we  define  the  potential 
energy  of  a  rigid  uniform  ring  of  radius  R  with  its  center  displaced  by  Y  from 
the  origin,  as 


VR’  Y»  «>  -  Pb(R)  If  d5i  lr  +  Y|  -  R)  ♦1(?,  a),. 


and  the  potential  energy  of  the  entire  beam  slice,  also  displaced  by  Y,  as 


*  (Y,a)  -  dR  2ttR  $  (R,Y,a). 
s  ;  o  r 


To  simplify  notation,  we  assume  that  Y  >  0.  We  now  proceed  through  a  series 
of  physical  arguments  that  specify  the  sources  and  the  sign  of  the  force  on 
beam  elements,  rings,  and  on  the  beam  slice.  Although  we  have  used  an 
expansion  in  small  Y  to  derive  the  separation  of  the  plasma  charge  into  p 
and  ?2»  Eq.  (11),  we  do  not  need  to  invoke  the  smallness  of  Y  in  some  of  these 
arguments . 

V  V  V 

(i)  The  part  of  p^(r)  with  r  >  r^  exerts  no  net  force  on  a  beam  element 
at  location  r^.  This  is  apparent  from  Eq.  (14). 


.  . . . .  ■ 


*.  --  •-  . 


(ii)  It  follows  that  a  ring  of  plasma  charge  with  radius  r  exerts  no  net 
force  on  a  ring  of  beam  charge  with  radius  R  if  r  >  R  +  Y,  since  the  beam  ring 
is  completely  contained  within  the  plasma  ring. 

(iii)  It  is  also  true  that  an  element  of  plasma  charge  exerts  no  net 
force  on  a  ring  of  beam  charge  with  radius  R  if  it  is  contained  within  the 
ring.  This  follows  because,  from  (i)  the  beam  ring  exerts  no  force  on  the 
plasma  element  and  according  to  Newton's  third  law,  the  reverse  is  true  as 


well. 


(iv)  It  follows  from  (ii)  and  (iii)  that  if  R  >  Y,  the  force  on  the  beam 


ring  of  radius  R  is  determined  only  by  the  plasma  charge  within  the  crescent 
shaped  region  between  r  *  R  and  r  =  R  +  y,  as  shown  in  Fig.  13.  (We  recall 
that  r  is  the  beam-centered  coordinate,  r  =  |r  -  Y|.)  It  is  obvious  then  that 
the  electron  beam  ring  is  attracted  back  toward  Y  =  0,  if  the  plasma  charge 
p^(r)  is  positive  throughout  R-Y<r<R+Y,  and  that  the  beam  ring  is 
pushed  toward  larger  displacement  Y  if  p^(r)  is  negative  throughout  • 
R-Y<r<R+Y,  i.e.,  that  3$r/3Y  is  positive  in  the  former  case,  and 
negative  in  the  latter.  If,  in  particular,  Y  is  infinitesimal,  then  the  sign 
of  the  force  3$r/3Y  is  determined  entirely  by  the  sign  of  p^(R).  When  we  note 
the  properties  (15a)  and  (15b)  of  p^(r),  it  is  apparent  that  3  <t  r/3  Y  is 
positive  for  R  <  rQ  and  negative  for  R  >  rQ,  i.e.  the  inner  portions  of  the 
beam  are  attracted  toward  the  channel  and  the  outer  portions  repelled,  if  the 
displacement  Y  is  small^®. 

Next  we  consider  the  force  34>g/3Y  on  the  entire  beam  slice,  for  the  case 
of  a  small  displacement  Y  of  the  beam  from  the  channel.  According  to  (22)  and 
(23), 


$  (Y,a)  -  /  dR  R  p  (R)  /  d8  «.(r  -  R  +  Y  sin  8), 

S  D  1 


r.i' 


Expanding  to  order  Yz  and  using  Eq.  (14)  this  becomes 


$  (Y,  a)  *  $  (0,  a)  -  Y2  2tt2  /°°  dR  R  p.  (r 

S  S  '  O  D 


R)  P1  (r 


R). 


(25) 


The  beam  slice  is  attracted  toward  the  channel  if  the  integral  in  (25)  is 

negative,  or  repelled  if  it  is  positive. 

If  the  slice  is  narrow,  i.e.,  a  «  1  and  falls  off  much  faster  than 

p^,  the  integral  is  determined  entirely  by  p^(r  <  r  )  and,  according  to  (21a), 

is  negative.  The  slice  is  attracted  toward  the  channel,  since  it  is  made  up 

entirely  of  narrow  rings  which,  we  have  seen,  are  attracted  toward  the 

channel.  In  the  opposite  limit  of  a  broad  beam,  the  p^(R)  can  be  approximated 

by  P^(0)  in  the  integral,  and  then  the  integral  is  positive,  according  to 

(21c).  In  this  case  the  central  part  of  the  beam  is  attracted  toward  the 

channel  and  the  outer  part  repelled,  but  the  latter  dominates.  The  transition 

2  2 

between  a  narrow  beam  with  3  $  (Y  *  0,  a)/3Y  >0  and  a  broad  beam 

s 

2  2 

with  3  <&  (Y  ■  0,  a)/3Y  <  0  occurs  for  the  value  a  of  a  defined  by 


/  dR  R  pb(R)  Pl(R)  -  0. 


(26) 


Using  Eqs.  (14)  and  (18)  for  p^  and  p^,  we  find  after  some  algebra  that  aQ  is 
defined  as  the  solution  of  the  equation 


(1  +  oQ)  (2  +  oQ)  -  2  ln[oo(2  +  aQ)(l  +  ao)“2]. 


(27) 


Next  we  consider  the  limit  of  a  beam  that  is  wholly  outside  the  channel. 


i.e.  Y  »  rb  and  Y  »  o' 


-1/2. 


In  this  case,  it  is  clear  from  two-dimensional 
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electrostatics  that  the  beam  will  polarize  the  channel,  leading  to  an 
attractive  force  between  beam  and  channel  that  is  proportional  to  Y  .  This 
holds  whether  the  beam  is  narrow  or  broad  compared  to  the  conductivity 
channel.  We  have  already  seen  that  in  the  other  limit  Y  0,  the  beam  channel 
force  is  attractive  for  a  narrow  beam  a  «  1)  and  repulsive  for  a  broad 
beam  (a  »  1).  Thus,  we  can  sketch  the  potential  $s  (Y;  a).  This  is  done  in 
Fig.  14.  If  the  beam  is  narrow  (Fig.  14a),  the  equilibrium  of  the  beam  is  on- 
axis  in  the  channel,  but  if  the  beam  is  broad  (14b),  the  equilibrium  position 
of  the  beam  is  displaced  from  the  channel  axis  by  some  quantity  YQ,  where 
YQ/r^  is  a  function  of  a. 
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IV.  Conclusions 

In  this  report,  we  have  studied  in  some  depth  the  conditions  in  which  a 
beam  can  track  a  preformed  channel.  The  results  can  be  summarized  as  follows: 

(a)  As  long  as  the  initial  channel  conductivity  a  (“  4-wor^/c)  is  less 
than  0.8,  the  beam  always  tracks  the  channel,  in  the  sense  that  the 
head  of  the  beam  is  attracted  to  an  equilibrium  position  that  is 
either  on  the  channel  axis  or  off-axis  in  the  channel  by  a  specified 
displacement  that  is  typically  of  tho  order  of  the  beam  radius. 

(b)  For  a  pure  density  channel,  the  equilibrium  is  always  off-axis. 

For  other  types  of  channels,  the  equilibrium  position  of  each  beam 
slice  is  at  the  channel  axis  or  off-axis  depends  on  the  ratio  of 
channel  radius  to  beam  radius,  r^/rc-  For  a  pencil  beam,  the 
equilibrium  is  off  axis  if  r^/rc  >  1.0.  This  may  mean  that  the 
equilibrium  is  typically  off  axis  for  a  self-consistent  modeled  beam 
with  an  expanded  head.  If  rb/rc  <  ^,0  for  a  pencil  beam,  the 
equilibrium  is  on  axis. 

(c)  A  pure  conductivity  channel  (cjq  -  0.1)  gives  the  best  on-axis 
tracking.  A  combination  of  conductivity  and  density  channel  gives 
the  strongest  off-axis  tracking. 

(d)  The  central  part  of  the  beam  always  feels  a  net  force  attracting  it 
to  the  channel,  even  when  the  average  force  on  the  whole  beam  slice 
is  detracking.  But  for  this  case,  the  portion  that  feels  a  tracking 
force  is  usually  a  small  part  of  the  beam  slice.  Since  the  tracking 
force  is  usually  at  least  10  times  smaller  than  the  local  pinch 
force,  which  holds  the  slice  together,  we  do  not  expect  the  beam  to 
tear  itself  apart  with  the  central  part  tracking  the  channel  and  the 
outer  part  expelled  from  the  channel.  Our  model  is  not  yet  adequate 
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to  test  this  supposition,  however.  (Nor  are  linearized  particle 
simulation  codes  adequate  for  this  purpose.) 

(e)  The  ratio  of  the  tracking  force  to  pinch  force  at  full  current 
(Ft/Fp)  is  about  10“^.  The  angle  of  deflection  due  to  channel 
tracking  is  therefore  very  small.  We  can  estimate  its  magnitude 
as 

9  -  (I/IAY)1/2(Ft/Fp)1/2.  (28) 

For  typical  laboratory  experiments,  this  tracking  force  is  too  small  to  be 
detected.  For  example,  an  NRL  experiment  was  performed  with  an  8  kA,  1  MeV 
beam  with  radius  0.75  cm,  propagated  into  a  reduced  density  channel  with 

*  0.25  atm  on  axis  and  very  little  residual  conductivity^.  No  appreciable 
movement  of  the  beam  centroid  was  detected  after  one  meter  of  propagation,  in 

_3 

agreement  with  Eq.  (28)  which  gives  a  deflection  angle  9  ~  1.5x10 
radians,  too  small  to  be  measured.  In  a  more  recent  experiment,  the  same  beam 
was  propagated  in  a  reduced  density  ammonia  channel  with  on-axis  density  10 
torr  and  background  density  40  torr.  It  was  found  that  the  beam  body  hosed 
violently,  and  no  tracking  was  detected. 
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Fig.  2  Typical  c -dependences  of  beam  current  density  on  the  beam  axis 


plasma  current  density  Jp,  axial  electric  field  Ez  and 


conductivity  a  on  the  channel  axis  near  the  beam  head. 
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Fig.  3  The  variation  of  the  tracking  force  F t(?)  vs  5  for  a  wide  chann 
(rc  ■  2.0)  and  a  small  beam  displacement  (Y  ■  0.1). 
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The  dependence  of  the  tracking  force  Ft(?)  on  various  beam 
displacements  Y  for  a  narrow  channel  (rc  *  0.5).  Curves  (; 
for  a  +  n  channel.  Curve  (a)  is  for  Y  ■  0.1,  curve  (b)  is 
0.5,  curve  (c)  is  for  Y  ■  1.0  and  curve  (d)  is  for  Y  *  1.5 
(e)  is  for  the  n  channel  with  Y  ■  1.5  and  curve  (f)  is  for 


-d)  are 
for  Y  - 
Curve 


the  a  channel  with  Y  ■  1.5 


TRACK 


x. 


The  dependence  of  the  tracking  force  as  a  function  of  beam 
displacement  Y  for  r_  ■  0.5  and  c  »  9.0.  Channel  type  is  a  +  n 


The  tracking  forces  for  two  values  of  channel  radius  rc  -  2.0  and  r 
■  5.0  are  compared.  The  beam  displacement  Y  *  0.1. 
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Fig.  11  The  dependence  of  the  tracking  force  on  the  on-axis  channel 
with  Y  *  1.5  for  a  density  channel.  Solid  curve  denotes  p^ 


density 
*  0.1  atm 


and  dashed  curve  denotes  p  *0.2  atm. 
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Fig.  12  The  tracking  force  ft  on  a  ring  of  beam  electrons  located  between  r 
and  r  +  ir  as  a  function  of  r.  Two  values  of  beam  displacements  are 
plotted,  Y  *  0.1  and  0.5.  The  central  part  of  the  beam  always 
tracks. 
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